Halperin-Saslow modes as the origin of the low temperature anomaly in NiGa 2 S4 
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The absence of magnetic long range order in the triangular lattice spin-1 antiferromagnet NiGa2S4 
[l[ has prompted the search for a novel quantum ground state. In particular, several experiments 
suggest the presence of a linearly dispersing mode despite no long-range magnetic order [H, 0, HI, 0] ■ 
We show that the anomalous low temperature properties of NiGa2S4 can naturally be explained 
by the formulation developed by Halperin and Saslow[f| where the linearly dispersing Halperin- 
Saslow mode may exist in the background of frozen spin moments and zero net magnetization. We 
provide highly non-trivial consistency checks on the existing experimental data and suggest future 
experiments that can further confirm the existence of the Halperin-Saslow mode. Our results place 
strong constraints on any microscopic theory of this material. 

PACS numbers: 



Introduction- Frustrated magnets are an excellent 
playground for the discovery of novel quantum phases 
of matter. The intricate interplay between geometrical 
frustration and quantum fluctuations lies at the heart of 
the physical mechanism for such possibilities. Recently, 
various frustrated magnets with small spin (S = | and 
5 = 1) have been discovered [H, @]. These experiments 
provide great opportunities for studying quantum as well 
as thermal fluctuation effects in frustrated magnets. In 
particular, NiGa2S4 is an insulator composed of spin-1 Ni 
atoms on a triangular lattice, that shows no long range 
magnetic order in neutron scattering down to at least 
0.35 K, despite an antiferromagnctic Curie-Weiss tem- 
perature |Ocw| = 80 K Instead, spin freezing is ob- 
served at Tf = 8 K, a temperature that also serves as the 
onset for a variety of unusual experimental signatures. 

Below Tf, the spin susceptibility saturates to a con- 
stant value, while the magnetic part of the specific heat 
acquires a T 2 power law [lj. This T 2 dependence of the 
specific heat is consistent with the presence of gapless hy- 
drodynamic modes propagating in two dimensions, pro- 
vided that these modes are coherent over a length scale 
L a that exceeds 500 lattice spacings On the other 
hand, elastic neutron scattering on powder [l[ and single 
crystal Q samples show no evidence of long range mag- 
netic order. Instead, coincident with the spin freezing 
temperature, there is the onset of static, yet short range 
magnetic correlations, with a correlation lenght £ of order 
of only 7 lattice spacings. In addition, nuclear magnetic 
resonance (NMR) measurements have observed a spin- 
lattice relaxation rate 1/Ti cx T 3 for T < Tf 0]. Recent 
inelastic neutron scattering measurements observe a lin- 
early dispersing excitation mode 0] 

In order to account for some of these results, there 
have been recent proposals of the formation of subtle 
types of long-range spin ordering. Based on models with 
large biquadratic interactions, Tsutentsugu and Arikawa 
proposed a three-sublattice nematic state 0] and Bhat- 



tacharje et al. proposed a uniform nematic state [9( (see 
also Ref.[ll| for related theoretical models). On the other 
hand, Kawamura and Yamamoto proposed a state with 
bound Z 2 vortices (lQj . The nematic states of Refs. [H 
and Q may be consistent with the constant spin suscep- 
tibility, and also contain gapless director modes, which 
can account for the T 2 specific heat. However, these pro- 
posals leave a few questions unanswered. For instance, 
there is no clear connection between the onset of nematic 
order and the appearance of static short-range magnetic 
correlations seen at Tf. In addition, a nematic state does 
not account for the linearly dispersing inelastic neutron 
scattering peak, or for the ~ T 3 spin-lattice re- 

laxation rate - the director modes of a nematic do not 
couple directly either to the neutrons or to the atomic 
spins probed by NMR[ll|. 



In this Letter, we will argue that these experimental 
discoveries can be explained by the presence of gapless 
Halperin-Saslow (HS) modes [5J. Halperin and Saslow 
originally proposed the existence of hydrodynamic modes 
in the context of spin-glasses. However, as discussed be- 
low, HS modes can exist in a much broader class of ma- 
terials. HS modes have a linear dispersion, and they 
couple directly to atomic spins and also to neutrons. 
Therefore, they can account for the full phenomenology 
of NiGa2S4 in a natural way. In the rest of the paper, 
we will provide various qualitative and semi-quantitative 
consistency checks, and also discuss future experimen- 
tal consequences of the HS scenario. We are not aware 
of any confirmed example of the Halperin-Saslow modes, 
therefore NiGa2S4 may present the first clear case for 
these long-sought low energy modes. Our discussion will 
be phenomcnological in nature, as it will not address the 
microscopic mechanism by which the system chooses the 
HS state, but will rather rely on general hydrodynamic 
considerations. This has the advantage that our results 
are not limited to NiGa2S4, but are easily adapted to 
other frustrated magnets that may be in the HS state. 
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Overview of the Halperin-Saslow theory - Halpcrin and 
Saslow showed that a system with static moments can 
support low energy hydrodynamic modes, even in the 
absence of periodic long range magnetic order [f| . Their 
formalism is very general and relies on two assumptions 
only: the presence of static moments in the (possibly 
metastable) ground state, and a finite spin stiffness to 
slow spatial deformations of the static spin texture. If 
the moments are non-collinear and have zero net magne- 
tization, the HS theory predicts the existence of three 
linearly-dispersing hydrodynamic modes. The clastic 
neutron scattering experiments in NiGa2S4 indicate that 
there are static helical spin textures at temperatures be- 
low Tf . The fact that these spin textures are short ranged 
is inconsequential to the HS scenario, provided that a fi- 
nite spin stiffness survives at long distances. 

For a system satisfying these assumptions, the free en- 
ergy of a weakly-perturbed state above the ground state 
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AF[m,e] = ^J2 [ ^K^+ftW 2 ]. (1) 

a = l J 

Here x 1S the spin susceptibility, p s is the spin stiffness, 
m a (r) is the local magnetization density, and 9 a are 
locally-defined rotation angles that describe the defor- 
mation in the excited-state spin configuration. Here m a 
and 9 a are canonically conjugate variables. In Eq. ((T|) 
we have assumed that the system is composed of weakly 
coupled layers of thickness d. 

Starting from Eq. ([T]) one obtains three degenerate po- 
larizations of spin waves with dispersion ujk = vk, with 
velocity 



(2) 



where 7 = g^s/ft- is the gyromagnetic ratio. In a more 
general case, x an d p s are tensors, leading to three non- 
degenerate spin waves, with polarization-dependent ve- 
locities Vj = l(Pj/Xj) 1/2 - 

Basic consideration: contrasting specific heat and elas- 
tic neutron scattering data - The presence of gaplcss HS 
modes is consistent with the T 2 specific heat in NiGa2S4. 
For linearly dispersing modes in two dimensions. 
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Elastic neutron scattering measurements on both pow- 
der and single crystal samples show that the freezing tem- 
perature Tf coincides with the onset of static short-range- 
ordered moments. Peaks appear at an incommensurate 
wave vector Q = (1/6 — e, 1/6 — e) . The peaks are static, 
but they are not wave vector resolution-limited, corre- 
sponding to short-range order over a length £ xy ~ 7a 
in-plane, and £ z ~ 6A out-of-planc. The wave vector Q 
is near the half-way point between the center and the cor- 
ners of the Brillouine zone, and the spin moments seem to 
be arranged in a helical pattern. This wave vector and 
spin structure would be expected from a system with 
a weak ferromagnetic nearest neighbor coupling and a 
large antifcrromangctic third-nearest neighbor coupling, 
J 3 3> — J\ > 0. Such a hierarchy of couplings is consis- 
tent with semi-empirical cluster calculations 12j , and a 
larg e J3 > is also supported by ab initio calculations 
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Note that the correlation length £ xy is much smaller 
than the scale Lq, i.e. the collective modes responsible 
for the T 2 specific heat exist over a much larger length 
scale than the helical regions seen in neutron scattering. 
Here we explore the possibility that below Tf, the system 
is in a quantum spin-glass-like state composed of helical 
domains that do not order at long distances, yet have 
a finite spin-stiffness. One possible picture is to coarse 
grain the system, and to think of each helical domain 
as some block degree of freedom. The residual interac- 
tion between blocks, of order Tf, is frustrated and disor- 
dered, leading to freezing behavior at Tf in the absence 
of long-range helical order. However, independently of 
the microscopic mechanism by which such a state may 
be formed, we focus on the generic properties that follow 
from the hydrodynamics of such a system. 

Two energy scales from the susceptibility and specific 
heat.- As a first semi-quantitative check, we would like 
to know whether the magnitude of the low tempera- 
ture susceptibility and specific heat anomaly are consis- 
tent with the HS scenario. In this scenario, these two 
quantities are related through the spin stiffness p s , since 
cm IT 2 oc 1/v 2 (x Xm/ Ps- From xm and cm/T 2 , we can 
extract two energy scales, 



(3) 



£1 



2g 2 p 2 B S(S + l) 

zxm/Na 
3n p C(3) k% 



fc B (113K) 



(4) 



Xm 



g 2 fx 2 B c M /T 2 



fc s (6.5K) (5) 



Here, £(3) = 1.202..., Na is the Avogadro's number, 
the sum runs over polarizations j, Vj is the collective 
mode velocity, and v = \^3a 2 d/2 is the volume of a unit 
cell where we have assumed decoupled layers of thickness 
d and a is the in-plane lattice spacing. Here, we have 
allowed for the possibility that the collective modes are 
coherent only up to some length scale Lq. Experimental 
fits yield a lower bound, Lq > 500a [lj]. 



where z 



6 is the coordination number of the lat- 



tice, and n p = 3 is the number of degenerate hydrody- 
namic modes. Above, we have used the experimentally- 
measured values, xm(T — * 0) = 0.0089 emu/mole and 
c M /T 2 = 2.6 x 10~ 2 J mole- 1 K" 3 0]. 

The energy E± is related to the effective interaction 
between a spin and its environment, whereas E2 corre- 
sponds to the spin stiffness extracted from measurements. 
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Note that, while E\ is comparable to the microscopic en- 
ergy scale |Ocw| = 80 K, E2 is much smaller. This 
is consistent with the HS scenario: the stiffness is ex- 
pected to be considerably renormalized relative to the 
microscopic energies, to a scale of order of the freezing 
temperature Tf = 8 K. 

We can provide a theoretical upper bound on the spin 
stiffness. This is done by computing the spin stiffness of 
a long-range-ordered helical state, ignoring the fact that 
spin correlations in NiGa2S4 are short ranged. Carrying 
out an analysis similar to Ref. [H, 
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where the sum runs over all vectors fL connecting lattice 
sites, fi x = fi-x for an arbitrary direction x, and V is the 
volume of the sample. Equation (|6|) is an upper bound 
on the stiffness of NiGa2S4, since it assumes an clastic 
deformation of the spin texture that is uniform across 
the sample. On the other hand, in the real material with 
disorder, it pays off to put the bulk of the deformation 
across weak links. Since the helical spin correlations in 
NiGa2S4 have wave vector near Q = (1/6,1/6), Eq. ([6|) 
gives the upper bound, 
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Using \Q C w\ = 80 K and the ratio Q J%/J 3 = -0.2, 
we obtain J3 w fcs(25K). Plugging in the experimen- 
tal values [l|, g = 2, \(S t )\ w 0.75, we obtain an up- 
per bound, p™£ x d = fcs(49K), which is indeed larger 
than £2- Conversely, we get a lower bound cm/T 2 > 
2.9 x 10~ 3 J mole _1 K~ 3 that is consistent with the exper- 
imentally measured value, which is about 9 times larger 
than this bound. In fact, the ratio between the two 
is comparable to the frustration parameter, Tf /\&cw\i 
which is a natural renormalization factor for the spin 
stiffness in a frustrated spin system. 

Note that the definition of the two energy scales E\ and 
E2 is model independent, and can therefore be used to 
characterize the experiments and to compare with other 
theoretical proposals. As we have argued, in the HS sce- 
nario, it is natural to find E2 ~ Tf -C E\ ~ |Qcw|- On 
the other hand, for the spin ncmatic states of Refs. Q 
and 0], this requires some fine tuning. For example, in 
order to obtain E2 -C |6cw|, the system has to be tuned 
to be very close to a quantum critical point. 

Effect of magnetic field and spin anisotropy.- One of 
the salient features of NiGa2S4 is the weak effect of mag- 
netic fields on the specific heat anomaly. For a typi- 
cal magnetically-ordered system, an applied field gaps 
out some of the collective modes and therefore sup- 
presses the specific heat below a temperature of order 
T* = gn B H/k B - 

In the HS scenario, an applied field splits the 3 polar- 
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Hence, ujq is unchanged by the magnetic field and lu + sa 
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g^BH + na ^ H is gapped. However, the third mode, 
obtains a quadratic spectrum at low en- 
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crgics. This soft mode compensates for the gapped mode 
w+, so that there is negligible deviation from the T 2 
anomaly down to T w T*/4 as can be seen in Fig[T] 
This is consistent with the weak magnetic field depen- 
dence of the specific heat down to the lowest temperature 
in the experiments. Below T ps T*/4, lu- dominates, 
leading to a low temperature regime with at T-lincar 
specific heat anomaly. Notice that nuclear quadrupo- 
lar resonance (NQR) at Ga sites shows the existence of 
inhomogeneous internal field. The resulting distribution 
of Zccman energy would further mask the magnetic field 
dependence of the specific heat. 

In realistic systems, anisotropy in the spin interactions 
will cut off the T-linear specific heat. Susceptibility mea- 
surements suggest that NiGa2S4 has a weak easy plane 
anisotropy. In the absence of a magnetic field, this type 
of anisotropy introduces a gap in two of the three HS 
modes. This gap is of order of the anisotropy energy, 
and is likely very small in NiGa2S4. 

Nuclear Magnetic Resonance - NMR measurements on 
69 ' 71 Ga nuclei yield a spin lattice relaxation 1/Ti oc T 3 at 
temperatures below 1 K 0. This is consistent with the 
HS scenario. The dissipative part of the dynamic spin 
susceptibility due to the HS modes is 
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plus regular terms. For modes propagating in two dimen- 
sions, this leads to 1/Ti oc T 3 . Note that, even if the HS 
modes arc gapped due to anisotropy, two magnon Raman 
scattering also yields 1/Ti oc T 3 [15}. On the other hand, 
for a nematic state, the dynamic structure factor acquires 
one more power of |k| at long wave lengths. Thus, for a 
nematic state we expect 1/Tj to vanish at least as T 4 at 
low temperatures, in contrast to experiment. 

Direct detection of HS modes via inelastic neutron scat- 
tering and future experiments - The HS modes are spin- 
one excitations, and therefore couple directly to neutrons. 
Recent inelastic neutron scattering measurements indeed 
see evidence of a linearly-dispersing mode centered at 
wave vector Q, and with spin- wave velocity f nC utron = 29 
mcV AQ. On the other hand, the spin-wave velocity ex- 
tracted from the specific heat anomaly via EqJ3] is much 
smaller, 

v s h = V^5.4 meVA = 9.35 meVA, (8) 
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T/T* 

FIG. 1: Specific heat for an isotropic system (thick solid line); 
for an isotropic system in a magnetic field (dashed line); and 
for a system with easy-plane anisotropy with anisotropy en- 
ergy 0.3T*, together with an magnetic field perpendicular to 
the anisotropy plane (thin solid line). Here, T* = gfisH . 

where the last value corresponds to n p = 3 collective 
mode polarizations. 

This large difference between v s h and fneutron is due 
to the fact that ^neutron probes the velocity of relatively 
high-energy spin-waves. In the current measurements, 
the inelastic neutron peaks can only be resolved clearly 
for energies larger than 1 meV. These relatively high- 
energy modes do not contribute directly to the specific 
heat anomaly below Tf = 8K. In fact, at these energies 
the spin waves are expected to probe the bare stiffness. 
An estimate for the bare stiffness tensor is the stiffness 
computed for the helical state, as given in Eq. ([7]). This 
yields two different spin-wave velocities, 

fz.holical = «j/,helical = 24 me V A (9) 

v z , helical = 34 me VA. (10) 

These velocities are consistent with current experiments 
- note that the peaks measured with neutrons are very 
broad, and thus it is very difficult to identify two separate 
but nearby peaks. 

Future experiments may be able to distinguish the two 
separate spin-wave velocities. In addition, with increased 
resolution it may be possible to measure the dispersion 
of spin- waves to lower energies. In the HS scenario, this 
would yield a spin-wave velocity that changes with wave 
vector, from v s -h very close to Q, to Uheiicai at wave vec- 
tors far from Q. Hence, this would constitute a direct 
measurement of the length-scale dependent spin stiffness. 
One further signature of the HS modes that may be seen 
in future neutron scattering experiments is the spin- wave 
damping, which is predicted to be [j| 

A = L>(k-Q) 2 (11) 

We note that the low energydirector modes that occur in 
the nematic states of Refs. do not couple directly to 



neutrons at long wave lengths [ll| , and therefore cannot 
account for these inelastic neutron scattering peaks. 

Other considerations: two level systems - In conven- 
tional spin-glasses, the T 2 contribution to specific heat 
coming from Halperin-Saslow modes is overwhelmed by 
a linear T contribution coming from localized two-level- 
systems (TLS) [3, 

ctis = ^-N a k 2 B T. (12) 
o 

Here, Nq is the density of states of TLS. A rough estimate 
for Nq can be obtained by assuming that the energy dis- 
tribution function of a single TLS is 1/ {ksTf), and that 
there is one TLS in each correlation volume v% ~ (,xyd- 
Using the in-plane correlation length £ xy ~ 7a, one ob- 
tains an approximate upper limit Nq < l/^ 2 y ^ 2 fcsTj, 
leading to 

< 3.0X10- 2 J mol^K- 2 . (13) 

If present, such a contribution would be seen only at the 
lowest temperatures explored in Ref. [l[. This suggests 
that the HS modes give the dominant contributions in 
the temperature range studied in Ref. [lj]. 

Conclusion - We have shown that the phenomenology 
of NiGa2S4 below the freezing temperature Tf can nat- 
urally be explained by the Halperin-Saslow scenario. In 
this theory, linearly dispersing HS modes in the back- 
ground of frozen spin moments can exist for a much 
longer length scale than that of the short-range magnetic 
order. This explains the simultaneous onset of the short- 
range helical magnetic order and T 2 specific heat at Tf. 
The HS modes are consistent with T 3 spin-lattice relax- 
ation rate observed in the NMR experiments. We also 
showed that the specific heat anomaly in the HS sce- 
nario is only weakly sensitive to applied magnetic field 
as observed in the experiments. In more recent experi- 
ments on purer samples [13] , indications of freezing start 
to show up at 10K. However, only at 2 K do the mag- 
netic moments seem to be completely frozen fl7j ]. In the 
temperature range between 2 and 10 K, the moments 
may not be completely static. Given that the freezing 
temperature may depend on the time scale of the probe, 
there might be some ambiguity associated with it, lead- 
ing to some uncertainty of the value of the spin stiffness 
in our analysis. Notice, however, that the value of the 
experimentally-extracted spin stiffness Ei = 6.5 K is in 
rough agreement with both extremes of the temperature 
range 2-10 K. 

It was pointed out that the HS modes may have already 
been observed in inelastic neutron scattering at relatively 
high energies, although current experimental resolution 
does not give access to the low energy mode dispersion 
relation. We argued that the spin-wave velocity mea- 
sured at relatively high energies reflects the bare spin 
stiffness and the velocity of the low energy HS modes, 
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that is directly related to the low temperature specific 
heat, should be different at lower energy scales. We pro- 
vided self-consistent estimates of both the low and high 
energy spin wave velocities. These predictions as well 
as the three polarizations of the spin waves and the de- 
tailed magnetic field dependence of the HS modes can be 
tested in future neutron scattering experiments on single 
crystals. 
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